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We present a new analysis of the helicity parton distributions of the nucleon. The analysis takes 
into account the available data from inclusive and semi-inclusive polarized deep inelastic scattering, 
as well as from polarized pp scattering at RHIC. For the first time, all theoretical calculations are 
performed fully at next-to-leading order (NLO) of perturbative QCD, using a method that allows 
to incorporate the NLO corrections in a very fast and efficient way in the analysis. We find evidence 
for a rather small gluon polarization in the nucleon, over a limited region of momentum fraction, 
and for interesting flavor patterns in the polarized sea. 

PACS numbers: 13.88.-|-e, 12.38.Bx, 13.60.Hb, 13.85.Ni 
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Introduction. — The exploration of the inner structure 
of the nucleon is of fundamental importance in Nuclear 
and Particle Physics. Of particular interest is the spin 
structure of the nucleon, which addresses questions such 
as hovir the nucleon spin is composed of the spins and 
orbital angular momenta of the quark and gluons inside 
the nucleon. In deep inelastic scattering (DIS) of Icptons 
off polarized nucleons it was found that surprisingly little 
of the proton spin is carried by the quark and antiquark 
spins pj. This has triggered much theoretical progress, 
and led to new experiments dedicated to unraveling the 
proton spin structure. Among them, experiments in po- 
larized proton-proton collisions at the BNL Relativistic 
Heavy Ion Collider, RHIC, have recently opened a new 
stage in this quest. 

The structure of a nucleon in a helicity eigenstate is 
foremost described by the (anti)quark and gluon helicity 
parton distribution functions (PDFs), defined by 



Af,ix,Q^) = f+ix,Q')-f-ix,Q^). 



(1) 



Here, fj'ix, Q^) [Jj{x, Q^)] denotes the distribution of a 
parton of type j with positive [negative] helicity in a nu- 
cleon with positive helicity, having light-cone momentum 
fraction x of the nucleon momentum and being probed at 
a hard scale Q. The integral A/j (g^) = /^ A/j (x, Q'^)dx 
measures the spin contribution of parton j to the proton 
spin, which is one reason why there are world-wide efforts 
to extract the Afj{x, Q^) from experimental data. 

The non-perturbative but universal A/j are accessible 
in measurements of double-spin asymmetries. 



A 



dAa da^ 



da+- 



LL 



da 



<ic^^ 



da+- 



(2) 



for processes characterized by large momentum transfer 
and helicity settings ±. Taking high transverse momen- 



tum {pt) reactions in polarized pp scattering as an exam- 
ple, the cross section at hadron-level schematically reads 
up to corrections suppressed by inverse powers oi px'- 

dAa ^ J2 h^^'' dxtAfa{xa,Q'^)Afb{xb,Q^) 



X dAaab{Xa,Xb,PT,a.s{Q'^),PT/Q)- 



(3) 



The sum runs over all initial partons a, 6, with dAuab the 
corresponding partonic cross sections, defined in analogy 
with Eq. ([2]). We note that depending on the experimen- 
tal observable, also an additional fragmentation function 
may occur in Eq. ([3]). An equation similar to ^ holds 
for the unpolarized cross section da. The dA&ab depend 
only on scales of the order of the hard scale px and are 
hence amenable to QCD perturbation theory. A consis- 
tent NLO analysis of ^ requires use of NLO partonic 
cross sections and scale evolution for the PDFs. 

In DIS, an expression analogous to the one in Eq. ([3]) 
holds, except that there is only one parton distribution. 
Efforts over the past three decades have produced ex- 
tensive data sets for polarized DIS [2|- Results from 
semi-inclusive DIS (SIDIS) [1 Q, IN -^ IhX, with h 
an identified hadron in the final state, have the promise 
to put individual constraints on the various quark fiavor 
distributions in the nucleon. Recently, the first precise 
(in part still preliminary) All measurements from RHIC 
have emerged [J), which are expected to put significant 
constraints on the helicity gluon distribution, Ag{xjCp), 
along with results from lepton nucleon scattering [5| . 

This paper presents the first "global" NLO analysis 
of the data from DIS, SIDIS, and RHIC in terms of the 
helicity PDFs. While there have been quite a few NLO 
analyses of the polarized DIS data in the past (see, e.g., 
Refs. [a, lZI, |8| ) , some of which [81] also take into account in- 



the remaining inverse Mellin integrals in ([5]) can be per- 
formed extremely fast. In fact, the calculation of dAcr 
becomes well over two orders of magnitude faster than 
for the "direct" method of computing it via Eq. ([3|) . This 
brings one to the kinds of computational speeds needed 
for a full inclusion of pp scattering data in a NLO fitting 
analysis without having to resort to approximations. We 
note that the calculation of the arrays just mentioned 
is a major computational challenge; however, we found 
that this procedure can be made very efficient and fast 
by using adaptive Monte-Carlo sampling techniques. 

Another crucial issue to be addressed in a global anal- 
ysis is the estimate of the uncertainties in the extraction 
of the various A/j , associated with either experimental 
or theoretical uncertainties. We pursue here an approach 
based on the use of "Lagrange multipliers" (LM) 8l.lll|. 
It has the advantage that no assumptions regarding, e.g., 
a quadratic behavior of the x^ function around the min- 
imum need to be made. Rather, one investigates how ^ 
varies as a function of a particular observable or variable 
xAfj{x, Ql) = NjX°'' (1 - x)'^' (1 + 7j\/5 + r]jx), (4) of interest. The uncertainty range is then defined by the 



formation from SIDIS, the full inclusion of RHIC data in 
the NLO analysis is a new feature. The fact that the cross 
section in the pp case, Eq. ([3]), is bilinear in the PDFs, 
its more complicated kinematic structure, and the overall 
high complexity of NLO partonic cross sections, present 
significant technical challenges in this endeavor. Based 
on a technique presented in [9|, we have now developed 
and commissioned a systematic procedure for performing 
full NLO analyses also when using pp scattering data. 

Global analysis and Mellin technique. — The idea be- 
hind a global analysis is to extract the universal PDFs 
entering factorized cross sections such as Eq. (j3|) by opti- 
mizing the agreement between the measured spin asym- 
metries for DIS, SIDIS, and pp scattering, relative to the 
accuracy of the data, and corresponding theoretical cal- 
culations, through variation of the shapes of the polarized 
PDFs. To be specific, we choose an initial scale for the 
evolution of Qo = 1 GeV and assume the helicity PDFs 
to have the following flexible functional forms: 



with free parameters Nj^aj, (3j,"fj,r]j. The strategy is, 
then, to evolve the distributions at NLO to the scales 
relevant to the various data points, to use the evolved 
distributions to calculate the NLO theoretical spin asym- 
metries at the kinematics of each data point, and to con- 
struct a x^ function testing the goodness of fit. The 
optimal parameters are then found by minimizing x^- 

The main technical challenge in our analysis lies in the 
inclusion of the RHIC All data. The x^ minimization 
procedure easily requires of the order of 10^ or more eval- 
uations of the rather complex NLO pp cross sections [10| 
for each data point, so that the computer time needed for 
a minimization directly on the basis of Eq. ([3]) becomes 
excessive. In Rcf. |9|, we devised a method that allows 
to overcome this problem by working in Mellin-moment 
space. This is the approach that we will also pursue here. 

In short, it amounts to expressing the PDFs in Eq. ^ 
by their Mellin inverse transformations to find, 

dAa = -t^E/ ^^ f dmAf:{Q')Afr{Q') 

X dXa dXbX~"x^"^dAa-ab{Xa,Xb,...), (5) 



region for which the increase A^^ in x^ above its low- 
est value is tolerable. In this way, one finds the largest 
possible range for predictions of a certain physical quan- 
tity that is consistent with a given Ax^ . Use of the LM 
method requires performing a huge number of fits, for 
which the efhciency of our Mellin method is crucial. 

Results of global analysis. — We now present results for 
our NLO global analysis of DIS, SIDIS, and (in part pre- 
liminary) RHIC data. The data sets we take into account 
[2, y ell's listed in Tab. |T] together with their respective 
X^ values [12] ■ Notice that we use new sets of fragmenta- 
tion functions (FFs) from Ref. [13|, which are consistent 
with all relevant unpolarized SIDIS and RHIC data en- 
tering the spin asymmetries analyzed here. Uncertainty 
estimates of the FFs 13| are propagated in the compu- 
tation of All and included in x^ as a theoretical error. 
We use the strong coupling as and unpolarized PDFs of 



Ref. [ij]. Other sets [15| give very similar results. 



where Cn,7n denote appropriate integration contours in 
complex n, m-space, see [9| for details. Here, the crucial 
point is that the information on the PDFs, contained in 
the Mellin moments A/" and A/™, defined as usual by 
A/"(Q^) = /q dx x"-~^ Afj{x, Q^), has been separated 
from the numerically tedious and time consuming part 
involving the dAcTab in the second row of Eq. ([5]). As 
the latter do not depend on the PDFs, their values can 
be computed prior to the fitting procedure, on a suitable 
array of moments n and m. Once this has been done. 



Rather than imposing the standard SU(2) and SU(3) 
symmetry constraints on the first moments of the quark 
and antiquark distributions, we allow for deviations: 

AU-AV={F + D)[l + esv(2)\. (6) 

AU + AV ~ 2AS = {iF - D)[l + esv(i)\, (7) 

where AJ^= [A// + A/^(Q§), F + D^ 1.269 ±0.003, 
3F — D ^ 0.586 ± 0.031 ]2|, and £su(2,3) are free param- 
eters. In total we have fitted 26 parameters [l6|, setting 
lu.d.s,g = in Eq. djl). Positivity relative to the unpo- 
larized PDFs of Ref. [1J| is enforced at Qo- In Fig. [1] we 
compare the results of our fit using Q — px to RHIC data 
from polarized pp collisions at 200 GeV ^], included for 
the first time in a NLO global fit. The bands are ob- 
tained with the LM method applied to each data point 
and correspond to the maximum variations for All com- 



TABLE I: Data used in our analysis [3, U, the individual 
X^ values, and the total x^ of the fit. We employ cuts of 
Q,PT>1 GeV for the DIS, SIDIS, and RHIC high-pT data. 



experiment 



data 
type 



data points 
fitted 



X 



EMC, SMC 


DfS 


34 


25.7 


COMPASS 


DfS 


15 


8.1 


Ef42, Ef43, E154, Ef55 


DIS 


123 


109.9 


HERMES 


DfS 


39 


33.6 


HALL-A 


DIS 


3 


0.2 


CLAS 


DIS 


20 


8.5 


SMC 


SIDIS, ;i± 


48 


50.7 


HERMES 


SIDIS, h^ 


54 


38.8 




SIDIS, 7r± 


36 


43.4 




SIDIS, K^ 


27 


15.4 


COMPASS 


SIDIS, h^ 


24 


18.2 



PHENIX (in part prel.; 
PHENIX (prel.) 
STAR (in part prel.) 



200 GeV pp, n° 
62 GeV pp, 7r° 
200 GeV pp, jet 



20 
5 
19 



21.3 
3.1 

15.7 



TOTAL: 



467 



392.6 



puted with alternative fits consistent with an increase of 
Ax^ = 1 or Ax^/x^ = 2% in the total x^ oi the fit. 

Our newly obtained antiquark and gluon PDFs are 
shown in Fig. [2] and compared to previous analyses [6|,la]. 
For brevity, the total Au + Am and Ad + Ad densities 
are not shown as they are very close to those in all other 
fits [a, lZI, y. Here, the bands correspond to fits which 
maximize the variations of the truncated first moments. 



A/. 



l,[a;„ 



'(Q^) 



Af,ix,Q^)dx, (8) 



at Q2 = lOGeV^ and for [0.001 - 1]. As in Ref. [| they 
can be taken as faithful estimates of the typical uncer- 
tainties for the antiquark densities. For the elusive po- 
larized gluon distribution, however, we perform a more 
detailed estimate, now discriminating three regions in x: 
[0.001 - 0.05], [0.05 - 0.2] (roughly corresponding to the 
range probed by RHIC data), and [0.2—1.0]. Within each 
region, we scan again for alternative fits that maximize 
the variations of the truncated moments Ag^'I^"""^^"''"^!, 
sharing evenly to Ax^- In this way we can produce a 
larger variety of fits than for a single [0.001 — 1] moment, 
and, therefore, a more conservative estimate. Such a pro- 
cedure is not necessary for antiquarks whose a;-shape is 
already much better determined by DIS and SIDIS data. 
One can first of all see in Fig.[2]that Ag{x, Q^) comes out 
rather small, even when cornpared to fits with a "moder- 
ate" gluon polarization (y, |8|, with a possible node in 
the distribution. This is driven mainly by the RHIC 
data, which put a strong constraint on the size of A^ 
for 0.05 ^ a; < 0.2 but cannot determine its sign as they 
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FIG. 1: Comparison of RHIC data Hj and our fit. The shaded 
bands correspond to Ax = 1 and Ax /x = 2% (see text). 
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FIG. 2: Our polarized sea and gluon densities compared to 
previous fits [y, |g] . The shaded bands correspond to alterna- 
tive fits with Ax^ = 1 and Ax^/x^ = 2% (see text). 



mainly probe Ag squared. To explore this further. Fig. [3] 
shows the x^ profile and partial contributions Ax| of the 
individual data sets for variations of A^^'t""^^"-^! . A nice 
synergy of the different data sets is found. A small Ag at 
X ~ 0.2 is also consistent with data from lepton-nucleon 
scattering [5|, which still lack a proper NLO description. 
The small x region remains still largely unconstrained, 
making statements about Ag^ not yet possible. 

We also find that the SIDIS data give rise to a robust 
pattern for the sea polarizations, clearly deviating from 
SU(3) symmetry, which awaits further clarification from 




FIG. 3: The x^ profile (a) and partial contributions A^? (b) of 
the data sets for variations of Ag^'^°-°^-°'^^ at Q^ = lOGeV^ 



TABLE II; First moments Afj' 



at Q^ = 10 GeV^ 



best fit 



Ax'^1 



: 0.001 

AxVx^ - 2% 



Aw + A« 


0.813 


0.793 


Ad + Ad 


-0.458 


-0.416 


Alt 


0.036 


0.028 


Ad 


-0.115 


-0.089 


As 


-0.057 


-0.006 


As 


-0.084 


0.013 


AE 


0.242 


0.366 



+0.011 
-0.012 
+0.011 
-0.009 
+0.021 
-0.020 
+0.029 
-0.029 
+0.010 
-0.012 
+0.106 
-0.120 
+0.015 
-0.018 



0.793 
-0.416 
0.028 
-0.089 
-0.006 
0.013 
0.366 



+0.028 
-0.034 
+0.035 
-0.025 
+0.059 
-0.059 
+0.090 
-0.080 
+0.028 
-0.031 
+0.702 
-0.314 
+0.042 
-0.062 



the upcoming W boson program at RHIC. A particularly 
interesting result emerges for the polarized strange quark 
distribution: a fit that excludes the SIDIS data prefers a 
negative As, but with SIDIS data included, As is forced 
to be positive for x > 0.02, in agreement with a recent 
LO analysis in [3|] . From the fit we find breaking param- 
eters esu(2,3) in (O, d?]) very close to zero, so that the 
first moment of As must be negative. Therefore, in the 
full fit. As turns negative at small x, gaining most of its 
area there. This is also visible in Tab. [ill where we show 
the full and truncated first moments of our PDFs. The 
behavior of As also leaves its imprint on the quark sin- 
glet, AE. Notice that below x ~ 0.001, where no data 
are available, the contributions to the full moments are 
determined by extrapolation of the distributions rather 
than constrained by the fit. 

Conclusions. — We have presented the first global NLO 
QCD analysis of DIS, SIDIS, and preliminary RHIC data 
in terms of the helicity parton distributions. A technique 
based on the use of Mellin moments of the PDFs allows 
to efficiently incorporate the data from pp scattering at 
RHIC in a full and consistent NLO analysis. We have 
found that the RHIC data set significant constraints on 
the gluon helicity distribution, providing evidence that 
A(7(x, Q^) is small in the accessible range of momentum 
fraction. We also found that the SIDIS data clearly point 



to a mostly positive Au and a negative Ad. The strange 
quark distribution As comes out negative at x < 0.02 
and positive at higher x, even though here the systematic 
uncertainties inherent in SIDIS are arguably largest. 

While our study should and will be improved on a 
number of aspects, in particular related to the inclusion 
of theoretical uncertainties and the treatment of experi- 
mental ones, we believe that it opens the door to finally 
obtaining a better and more reliable picture of the spin 
structure of the nucleon. In particular, it will help RHIC 
to realize its full potential, as hopefully more and more 
precise data will emerge over the next few years. We 
finally note that use of our fast and efficient Mellin tech- 
nique for incorporating NLO pp scattering cross sections 
in the analysis is of course not restricted to RHIC, but 
could equally find important applications at the LHC. 
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